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Decoupling is an important tool to prolong the coherence time of quantum systems. Most de¬ 
coupling schemes have been assuming selective controls on the system and it is believed that with 
global pulses one can only decouple systems with certain coupling terms like secular dipole-dipole 
coupling. In this article we show that with global pulses it is possible to reduce the coupling strength 
of other types of coupling, which we demonstrate with Ising coupling. The complexity of such pulses 
is independent of the size of system. 


I. INTRODUCTION 

Quantum systems suffer from decoherence due to interactions with environments. The task of decoupling is to 
remove unwanted couplings between systems and environments SQIil. Many decoupling schemes have been 
proposed and demonstrated in experiments for example, randomized dTOamical decoupling [l9j | 

uses randomly selected pulses at regular intervals, UDD (Uhrig dynamical decoupling) Q can cancel dephasing 
of a single qubit up to order n by using a minimal number of n pulses, CDD (concatenated dynamical decoupling) 
constructs decouplin g pu lse sequences recursively Q. There are also studies on using pulses to remove internal couplings 
of quantum systems |2ll . [^ . [37j | or engineer Hamiltonians BUS- 

A common feature of these decoupling schemes is that they all assume selective controls on the system. For many 
quantum systems, selective addressing of each qubit could be very high demanding, especially of those systems whose 
environment consists of the same physical objects as the system,for example, in some solid state devices, the system 
and the environment can be the spin of same nuclear species. For such systems, selective controls on the system is 
very hard, as the pulses usually affect all the spins, i.e., the pulses will be global. The known examples of decoupling 
with global pulses are WAHUHA [s^, MREV-8 and MREV-16 in Nuclear Magnetic Resonance, which exploit 

the symmetry of homonuclear secular dipole-dipole coupling to decouple the system. Such decoupling schemes rely 
on the symmetry of secular dipole-dipole coupling and therefore do not apply to other types of coupling. Recently 
applications of global pulses in Hamiltonian engineering were also studied[l8|. In this paper, we examine the use 
of global pulses to decouple the system with Ising coupling between the qubits, and show that, to the contrary of 
previous belief it is possible to decouple system with couplings different from secular dipole-dipole coupling. The 
advantage of global pulses is that the number of pulses needed for decoupling will be independent of the number of 
qubits, i.e. the complexity of global pulses is 0(1). 


II. AVERAGE HAMILTONIAN 

The principle of decoupling can be illustrated by the average Hamiltonian theory[^|3^, i.e., the propagator can 
be written as a single exponential relying on some Average Hamiltonian H which has the same effect as a time 
varying Hamiltonian. The full advantage of this theory is often realized in an interaction frame of a period and cyclic 
Hamiltonian. Assume that in an appropriate interaction frame, the Hamiltonian is piecewise constant Hi, H 2 , ■ ■ ■ Hm 
in corresponding time intervals ti,t 2 , ■ ■ ■ ,tm, then 


g-iHt _ 




Here Hi ■ ■ ■ Hm are transformed Hamiltonians from the physical Hamiltonian by applying pulses on the system, 
i.e., Hi = UjHUi where Ui represents the propagator generated by the pulses. The first few orders of the average 
Hamiltonian are 


H = H 


( 0 ) 
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FIG. 1: Spin topology: (a)Linear Chain (b)Square Lattice 


^ — -{Hiti + H2t2 + • • • + 

+ [^^3^3, ^^ 2 ^ 2 ] + •••}, 

^ = ^^{[^2^2, [^^2^2, ^^l^l]] — [^^2^2, ^^l^l]] + ■ • • } 

where t = 

A global pulse transforms an initial Hamiltonian H to 

Hi = {U ®U ®---® U)'^HU ®U ® ■ ■ ■ 

where U G SU (2) represents the propagator generated by the global pulse on each qubit. Global pulses have been used 
to decouple systems with secular dipole-dipole coupling. In this case, the coupling Hamiltonian takes the following 
form, 


Hf = ^ - IjJk. - Ijyhy), (2) 

jk 

here Ix ■= 5(10), ly '■= M 0 -?i) > the Pauli spin matrices and we denote Ity as the operator 

that acts as on the fth spin(see [ 13 ). 

Applying global (x) ^ pulse and (y) ^ pulse on the system, one gets 


Hf = Y, d^\‘^Ijyhy - Ij.Ik. - IjJk.). 

jk 

Ht - Ijyhy). 

jk 


( 3 ) 


It is easy to see that averaging these three Hamiltonians decouples the dipole-dipole coupling, i.e., 

St St St -mtH 

e^ezgj=g ^ 


where H = 0. This is the basic building block of WAHUHA, MREV -8 and MREV-16 and the effectiveness of 

However such decoupling scheme only works for 


such decoupling scheme has been experimentally demonstrated 
the systems with coupling al^Ix + Plyly + 744 where a -I- /3 -I- 7 = 0. 


III. REDUCE ISING COUPLING WITH GLOBAL PULSES 

At a first look, it may seem impossible to decouple Ising coupling with global pulses, as global pulses can not change 
the signs of Ising coupling. We will show that it is indeed possible by extending our previous study on Homonuclear 
decoupling [s^- 

Consider a system consisting of N qubits connected by Ising coupling, the coupling topology can take various 
shapes, for example, it can be a spin chain or spin lattice, as shown in Eig. O- 

A gradient magnetic field is added upon the system which induces Zeeman splitting on the qubits. The magnetic 
field and its gradient are large enough so that the differences of Zeeman splitting between coupled qubits are much 
larger than the coupling strength between the qubits. Under such gradient magnetic field, the Hamiltonian of the 
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system takes the form 


N 

^ — 'y ^ ^ y ^ ^jk^jz^kz (4) 

3 = 1 Uk)&G 

where ujj = -jl-Bj, fl is the magnetic moment and Bj indicates the magnetic field at site j. G is a graph indicating 
the coupling topology of the system, i.e., if the edge (jk) G G, then the qubits at site j and k are coupled. We assume 

iOi Jjk-i \^3 ^k\ Jjki'^i-i jtk. 

We will first use a two-qubit system to illustrate the decoupling strategy, then generalize it to N qubits with various 
coupling topologies. 

For a two-qubit system, the Hamiltonian is 

Hq = LUiIiz + UJ 2 I 2 Z + Jhzhz- (5) 

The basic pulse sequence consists of four periods, within each period it evolves according to the following Hamilto¬ 
nians 


Hi = COiIiz + OJ2I2Z + Jllzhz + + l 2 x)^ 

H2 = —UJiIiz — UJ2I2Z + Jhzhz + ^{Ilx + l2x), 
H 3 = —LOiIiz — LO 2 I 2 Z + Jhzhz — A{Iix -\- hx), 
Hi = Wi/iz + UI2I2Z + JIlzl2z — -^illx + l2x), 


where Hi is obtained simply by applying a magnetic field in x direction with effective amplitude A, H 2 is obtained 
from Hi by conjugating a 7r-pulse along the x direction, i.e., H 2 = Hi and H^ are obtained 

with a control field along the —x direction and conjugation with 7r-pulses along x direction. The 7r-pulses here are 
assumed to be infinitely narrow pulses. 

Each of the four Hamiltonian is maintained for a period of At, keep the terms up to the second order of the average 
Hamiltonian theory, we obtain the following average Hamiltonian over an interval of 4At, 

AAt 

Hleff — JIlzl2z T ^ ^2^2y') AAtJ(^Iiyl2z T ^lz^2y) 

(AAfV 4 

H-^-(wi/l2 -I- ^2/22) + J{^ly^2y — hzl2z) + 0((^At)^) 

Denote 9 = AtA and choose At such that 6* <C 1, but 0|wi — W 2 I J- Next we apply tt pulses along y direction to 
flip the signs of the third and fourth terms in Hie//. As a consequence, we create H 2 eff such that 


g-*ff 2 e// 8 Ai _ g-ii?ie// 4 AtgZ 7 r(/i;,+/ 2 y)g-iffie// 4 Atg-i 7 r( 7 i,,-|-/ 2 a) 

It is straightforward to see that 

H 2 eff — + ^ 2 l 2 y) + JIzSz + J{hyl 2 y ~ hzhz) + 0 { 9 ^). (9) 

The pulse sequence is shown in Fig. [5J 

The newly created Zeeman-like terms along the y-direction are orthogonal to the Ising coupling terms. Since we 
have 9\u3i — a; 2 | J, we can use rotating wave approximation to reduce the effective Hamiltonian to 

H^eff — + W2/2y) + -^9^JIiyl2y + 0{9^) (10) 

that is, if we repeat the procedure k times such that 8kAt 3 > , then [g-»^ 2 e// 8 Atjfe _ g-j 7 / 3 e// 8 feAt_ Compare 

with the original Hamiltonian Hq, we have effectively reduced the coupling strength by a factor of If we can 
create |wi — W 2 I ~ 10^ J, then 9 can be taken around ~ in this case i.e., the coupling strength is reduced 

by « 300 times. Further reducing of the coupling strength can be achieved by more iterations of the above procedure. 
Note that the local terms, such as -l-a; 2 / 2 y) can be canceled by Hahn echo pulses, i.e., inserting 7r-pulse along 

the X direction. 
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FIG. 2: Basic pulse sequence for decoupling. 


This decoupling strategy can be generalized to N qubits with various coupling topologies, which can be easily seen 
by substituting the Hamiltonian for A^-qubit system in Eq.(|6]), 


N N 

Hi = Jjkijzikz + Aljx, 

i=i (ifc)eG j 


N N 

H2 = — ''^^UJjljz + y]] Jjkijzikz + Aljx, 
1=1 (ifc)eG j 

N N 

II3 — ~ + y]] Jjkijzikz — y]] I^Ijx, 

i=l Uk)€G j 

N N 

H4 = y^ ^jijz + y^ Jjkijzikz— y^ Aijx, 


(11) 


i=l Uk)&G j 

Again by inserting 7r-pulses along ay direction and using rotating wave approximation, this creates an effective 
Hamiltonian 

n N ^ 

Heff = -j^'^^jljy + 2 ^^ y]] Jjkijylky + +0{9^). (12) 

i=l 0fe)6G 


The subsequent steps are similar as the ones outlined in the case of a two-qubit system. 

Since at every step of the analysis, the precise knowledge of the coupling strength is not required as long as they 
are small compare to the Zeeman splitting, this decoupling scheme is actually robust to the random fluctuation of 
the coupling strength Jjk- For example on the two-qubit system, the coupling strength during the four periods of the 
basic pulse sequence in Eq.(IS]) can be different from each other—we assume in each period the coupling strength is a 
constant drawn from a steady distribution, this assumption holds when the frequency of fluctuation is small compare 
to i—in this case, the basic pulse sequence becomes 


Hi = UJlIlz + UJ2I2Z + Jlllzl 2 z + A{Ilx + l 2 x), 
H2 = —Uliliz — UI2I2Z + J2I1ZI2Z + A(Jia; -b l2x), 
Hs = —Uliliz — OJ2I2Z + Jshzhz — A(/ia; + l 2 x), 
H4 = UJiIiz + UJ2I2Z + Jillzl2z — A(/ia; -b l2x)j 


(13) 


Keep the terms up to the second order of the average Hamiltonian theory, we obtain the following average Hamil- 
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tonian over an interval of 4At, 

>/i + ■/2 + <^3 + Ji 


H'uff = 


A/Si 

Ilzl2z H-+ UJ2l2y) 


+ 3 J 2 + 3i/3 + J 4 ) / ^ ^ T T \ (^At)^ 

H-g- {Ilyl2z + Ilzl2y) H-^+ t^ 2 ^ 2 z) 

^ ^ Ji + Ay, + J4 ^ + + 0((AA()») 


( 14 ) 


which reduces to Eq. © when the four coupling strengths are the same. Again by inserting tt pulses along y-direction, 
we can flip the signs of the third and fourth terms in and creates an effective Hamiltonian 


3-®^2e//8Ai _ g-iffh//4Atgi7r(/i^+/2a)g-iffh//4Atg-i7i-(/i^+/2^) 


(15) 


Here the two may have different coupling strength, we assume the coupling strengths for the 4 periods of 

basic pulse sequences for each are Ji, J 2 , ^ 3 , J 4 and J 5 , Je, J 7 , Js respectively. Then 


, Ji + J2 + ds + J4 + J5 + Je + T7 + ds . . 

'"26// — -^-tizt2z + 2 + ^2''2yj 


Ji + 3 J 2 + 3 J 3 + J 4 — J 5 — 3J6 — 3 J 7 — Js 


7Ji + J2 + J3 + 7J4 + 7J5 + J6 + J7 + 7J8^2 


24 


d{Ilyl2z + Ilzl2y) 

O^hyhy - hzhz) 


(16) 


0((AAt)3) 


where 6 * = AAt. Similar to the case with constant coupling strength, we choose At such that 0\uji — 022 ! ^ Ji- By 
rotating wave approximation the yz, zy and zz coupling are effectively averaged out as they do not commute with 
uJihy + ui 2 l 2 y, only yy coupling remains, whose strength is reduced by the order of 9^ comparing to the original 
coupling strength. The generation to N-qubit system is straightforward, similar to the case with constant coupling 
strength. 

This decoupling pulse sequence can also reduce the dephasing effect caused by the environment. Suppose each qubit 
in the system is coupled to environment where the coupling Hamiltonian is modeled as Hsb = f^<^z{gkbl. + g\bk)^ 
where b\,bk are bosonic operators for the fcth field mode of the environment, characterized by a generally complex 
coupling parameter The tt pulses along the x and y directions in our decoupling scheme also average out the 

net effect of Hsb- 

Finally, we present a numerical simulation illustrate the effects of our decoupling strategy. As shown in Fig.[3l the 
vertical axis represents the fidelity of U with respect to the identity operator, where fidelity measurement between 
two unitary operators Ui and U 2 is defined as 


\tr{U2Ul)\\triU,Ul)\y 

Other fidelity measures, for example, the average gate fidelity [H, can also be used, which is equivalent in 
our case. The simulation was done with 4 spins on a square lattice with one iteration of the pulses. Assume the 
lattice is put at a vicinity of a dysprosium micro-magnet with length of of 400/rm, width 4/tm and height 10/rm, 
which can generate a field gradient of dB/dz = 1.4T/rm“^[4l|. In addition a large homogeneous field Bq of ~ 7T 
is superposed. The distance of two neighbor nuclei spin is about ^ Inm, for the simulation we take the distance as 
Inm and the nuclear spin as ^sSigll. If the magnetic field gradient is put along the direction of y = ^x, then the 
Zeeman splittings for the 4 spins are 62.8kHz, 95.9kHz, 120.1kHz, 153.18kHz respectively. The secular component of 
the dipolar Hamiltonian which couples the ith spin to the /th spin is written 


= 


Att 


,1 — 3 cos^ 9i 


■BJi = 


where rp is the length of the vector connecting the spins and Op is its angle with the applied field. With chosen 
parameters, the coupling strengths between adjacent spins are ~ 17.3HZ, and ^ 6.1Hz between diagonal spins. Due 
to vibration of atoms, the couplings strength can be fluctuating, so we assume the coupling strengths at each instant 
of time follow independent normal distributions with mean values J = [17.3,17.9,18.5,19.2,6.1, 6 . 6 ]Hz and variance 
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FIG. 3: (Color online) The simulation was on a square lattice with four qubits. The coupling strengths between these 
qubits are fluctuating at each instant of time: the coupling strengths between neighboring qubits have mean equal to J = 
[17.3,17.9, 18.5,19.2]Hz, and variance equal to the J — 9 Hz; the coupling strength between diagonal qubits have mean equal 
to J = [6.1, 6.6]Hz, and variance equal to 3Hz. The Zeeman splitting of the four qubits are [62.8, 95.9,120.1, 153.18]kHz 
respectively. At = 10“^, 8 — A ~ 27r8 x 10® = 8kHz. The initial state of the qubit is . The blue line shows the 

fidelity without the global pulses, the red circles show the hdelity with the application of the decoupling pulses. 


equal to 10 Hz for the couplings between adjacent spins and 5Hz for the coupling between diagonal spins. At is taken 
to be 10“^s and A Ri 8kHz, 9 = ^. It can be seen that the global pulses reduce the decoherence rate by about two 
orders of magnitude in one iteration. Note that the results does not depend much on the precise numerical values, as 
long as the condition uJi ^ Jjk, \ujj — uJk\ ^ k is satisfied. 

IV. CONCLUSION 

We presented a method that reduces the Ising coupling strength of register qubits using global pulses. This can be 
used to reduce the residual coupling of quantum memories, where selective addressing may be hard or undesirable. 
The advantage of these global pulses is that the number of pulses does not grow with the number of qubits, i.e., the 
complexity of these pulses is only 0(1). This opens new directions for using global pulses for decoupling. 
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